Pump-probe coupling of matter wave packets to remote lattice states 
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The coherent manipulation of wave packets is an important tool in many areas of physics. We demonstrate 
the experimental realization of quasi-free wave packets of ultra-cold atoms bound by an external harmonic trap. 
The wave packets are produced by modulating the intensity of an optical lattice containing a Bose-Einstein 
condensate. The evolution of these wave packets is monitored in-situ and their reflection on a band gap is 
observed. In direct analogy with pump-probe spectroscopy, a probe pulse allows for the resonant de-excitation 
of the wave packet into localized lattice states at a long, controllable distance of more than 100 lattice sites from 
the main component. This coherent control mechanism for ultra-cold atoms thus enables controlled quantum 
state preparation, opening exciting perspectives for quantum metrology and simulation. 

PACS numbers: 03.75.Lm, 37.10.Jk, 67.85.-d 



Wave packets - non- stationary superpositions of energy 
states - have held a special place in physics ever since 
Schrodinger showed that a superposition of states can form 
a spatially localized entity evolving similar to a classical par- 
ticle [1]. Despite being ubiquitous in nature, wave packets 
remained an elusive concept for decades. The advent of ultra- 
short laser pulses has, however, enabled a multitude of ob- 
servations of wave packet phenomena in e.g. highly excited 
atomic Rydberg states 121 Ij], molecular states \^\m^ and semi- 
conductor systems f8f]. 

One of the main applications of wave-packets has been 
the investigation of complex quantum systems using so-called 
pump-probe spectroscopy y, [tI |8|]. An initial short pump 
pulse is used to create a wave packet in an excited state. Sub- 
sequently, the wave packet is allowed to evolve for a variable 
duration. Finally, a probe pulse interrogates this dynamic evo- 
lution by transferring it to a chosen final state. This process 
can thus provide information about coherent dynamics in the 
excited state Ig, Ml, but it can equally well be used for e.g. 
interferometry (lOj or for the preparation of a desired final 
state [11]. 

In parallel, the precise control of matter waves in Bose- 
Einstein condensates (BECs) has enabled a new field of wave 
packet manipulation [12]. For instance, Bragg diffraction [13y 
has been utilized to split and recombine wave packets interfer- 
ometrically [14] and freely propagating atomic wave-packets 
have been coupled out of a trapped BEC 1 15]. These imple- 
mentations, however, lack the structural complexity due to the 
combination of bound and continuum states present in e.g. 
molecular and semiconductor systems. 

In this letter we present the experimental realization of 
quasi-free wave packets of ultra-cold atoms created from an 
initial bound state. We produce these wave packets from an 
ultra-cold cloud by briefly modulating the intensity of an op- 
tical lattice in the presence of an external harmonic trap. The 
ensuing oscillatory motion of the wave packets is observed 
in-situ and shows only limited spreading. In an important ap- 
plication of this new phenomenon we couple the wave packets 
into localized lattice states |l6l llTO separated from the origi- 
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FIG. 1. (Color online) Spectrum of the ID single particle 
Schrodinger equation for a combined harmonic (cjtrap = Stt • 
37.9 Hz) and periodic potential at a lattice depth s = 16. The dashed 
red line indicates the maximum of the periodic potential. Using a 
short pulse of lattice amplitude modulation, quasi-continuum wave 
packets are created (1), which oscillate in the potential (2). An ap- 
propriately timed second modulation pulse can de-excite the wave 
packets into localized states (3). 



nal cloud by more than 100 lattice sites using a pump-probe 
approach directly analogous to the method discussed above. 

The production of wave packets and the population of out- 
lying localized states can best be visualized using the popula- 
tion density of the energy eigenstates. Figure [T] shows the cal- 
culated ID energy spectrum of the single particle Schrodinger 
equation in the combined harmonic and optical lattice po- 
tential d^ [nl]. It clearly shows pseudo-bands of allowed 
energies caused by the optical lattice, which are shifted ac- 
cording to the spatial dependence of the harmonic trap. Due 
to the site-to- site energy shift induced by the harmonic con- 
finement, tunneling is inhibited outside the trap center in the 
bound bands, resulting in so-called localized states. Above 
the lattice threshold (red dashed line) a quasi-continuum of 
delocalized eigenstates of the combined potential is formed. 
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FIG. 2. (Color online) Propagation of the wave packet after excita- 
tion with a 500 /is modulation pulse, a) Position of the upper ( x ) 
and the lower (o) wave packet as a function of the evolution time af- 
ter modulation, b) Wave packet dynamics obtained from a numerical 
simulation (see text). 



To realize this situation, the depth of the ID optical lat- 
tice ^Lat is chosen such that the second band is just barely 
trapped in the lattice, while higher bands are not. A harmonic 
modulation of the lattice depth with a frequency i^mod close 
to the transition frequency from the ground state to the sec- 
ond excited band, 1^02, can induce transitions between these 
bands [18] and further into the quasi-continuum states that ex- 
ist above the threshold of the lattice. 

In the following, we show experimentally that a short pulse 
of lattice modulation can induce this two-phonon transition 
into the quasi-continuum. This creates wave packets of nar- 
row quasi-momentum distribution which propagate outwards 
along the direction of the optical lattice. We also show 
that detailed knowledge of the dynamics can be used to in- 
duce a resonant de-excitation into the outlying localized states 
analogous to pump-probe spectroscopy. These processes are 
schematically indicated in Fig.[Tl 

Our experimental system [19] consists of a BEC in a har- 
monic magnetic potential and a one-dimensional optical lat- 
tice. A BEC of ^'^Rb atoms in the 5'^Si/2\F = 2,mF = 2) 
state is prepared in a quadrupole-Ioffe configuration (QUIC) 
trap [20] . After condensation, the magnetic potential is con- 
siderably relaxed to obtain trap frequencies of 12.3 Hz and 
37.9 Hz in the axial and radial direction respectively. Pure 
condensates of approximately 3x10^ atoms, are then loaded 
adiabatically into the vibrational ground state of a ID optical 
lattice. The optical lattice is formed along the vertical axis 
(one of the radial directions of the trap) by a retro-reflected 
laser beam at a wavelength of ALat = 914 nm with a 1/e^ 
waist of 120 fim. For all experiments the lattice depth was 
s = V^s^t/Erec = 16, where E^ec = h'^ /{2mXl^^). 

To produce the wave packets experimentally, a short 
(500 /is), large amplitude (0.3li,at) lattice modulation pulse 



is applied at a modulation frequency i^mod = 31 kHz close 
to z/02 = 30.5 kHz. Subsequently, images of the resulting 
clouds are taken in-situ following a variable evolution time. 
Two wave packets, traveling in each lattice direction away 
from the original cloud are observed as shown in Fig.[2]a) (in- 
set) [21]. The wave packets each contain roughly 2 % of the 
total atoms and do not show considerable spreading as they 
travel to the turning point, are reflected, and return. Due to 
gravity the dynamics of the wave packet traveling upwards is 
slightly modified compared to the one traveling downwards. 
For the quantitative analysis in the following we therefore use 
the upper component only. Figure [2] a) shows the distance of 
the wave packet from the center of the main component as a 
function of the evolution time during the first half oscillation 
period. It reveals a time to reach the turning point of roughly 
^turn = 2.6 ms (Tturn = 3.0 ms) for the upper (lower) peaks. 
We stress that this observation is clearly inconsistent with har- 
monic oscillation at the trap frequency entrap = 27r • 37.9 Hz 
from which we would expect a turning point at 6.6 ms. 

An intuitive understanding of the process can be obtained 
from an analysis based on Fig. [T] and the conventional band 
structure represented in the first Brillouin zone. The mod- 
ulation pulse transfers a small fraction of the ground state 
population via a two-phonon excitation to a superposition of 
quasi continuum states, determined by the individual coupling 
matrix elements. In the band picture, this is equivalent to 
a transfer into the fourth band for the chosen lattice depth. 
Energy conservation dictates a transfer into states with oppo- 
site quasi-momentum ifig'max, where 4 < qma^/h < 5 and 
kr = 27r/ALat- The wave packets exchange their initial ki- 
netic energy for potential energy by moving outward in the 
harmonic potential while rolling down the appropriate bands. 
As the first significant band gap is encountered, they inverse 
their direction due to a Bragg reflection. 

This interpretation directly allows for a simple calculation 
of the expected behavior of the wave packets. Treating the 
quasi-momentum hq as the momentum p of a classical particle 
moving in a harmonic trap (i.e. p{t) = hqmax cos(co'trap^)), 
the times at which the quasi-momentum passes the avoided 
crossings between bands n and (n — 1) at g = nkr are 

tn = arccos ( ^^^^ ) /co'trap for < n < 4. We obtain 

y <?max J ^ 

ts = 3.0 ms for qma^ ~ 4/cr in good agreement with the ex- 
perimental result and therefore interpret the reflection of the 
wave packet as a Bragg-reflection on the gap between the sec- 
ond and the third band (6-photon Raman process). A calcula- 
tion of the Landau Zener tunneling probability at the n = 2 
to n = 3 band gap supports this, since it is essentially zero at 
s = 16. 

For a more quantitative understanding, we have also per- 
formed a single-particle numerical simulation of the process 
by solving the time dependent Schrodinger equation using the 
rotating wave approximation. We do this by splitting the spec- 
trum into three parts: the zeroth band, the first and second 
excited band manifold, and finally the continuum states. Cou- 
pling between the individual states of separate manifolds is 




FIG. 3. Time dependence of the pump-probe spectroscopy, a) Pop- 
ulation in the locaUzed states as a function of the delay between the 
500 /is pump and probe pulses (z^mod = 32 kHz). The double peaked 
structure corresponds to the passage of the wave packet through the 
resonance region. Representative error bars are included at t = 1.75 
ms and at t = 6.75 ms. b) Numerical simulation of the pump-probe 
experiment (see text). 



introduced according to their coupling matrix elements. Al- 
lowing for this dynamics to evolve for the duration of the 
modulation pulse we obtain the excited wave packet as a su- 
perposition of eigenstates. The ensuing wave packet dynam- 
ics, dictated by the relative phase evolution of these states, 
is shown in Fig. [2]b) and reveals a turning point at 3.2 ms. 
Our simulation neglects interactions and since the mean field 
energy tends to flatten the effective potential we expect the in- 
clusion of interactions to decrease Tturn and thus give better 
agreement with the experiment. Furthermore, the simulation 
starts in the absolute single particle ground state, since inter- 
actions are neglected. We found that the simulation had to be 
performed at a relatively low modulation frequency (28 kHz) 
to obtain agreement with experiment and attribute this to the 
interaction induced broadening of the true initial state. 

The production process of wave packets described above 
bears strong resemblance to the pump pulse in a pump-probe 
experiment. We therefore illustrate the dynamical control of 
these wave packets by coupling them into outlying lattice 
states with a second, so-called probe, pulse. As mentioned 
above, these localized states appear since tunneling in the lat- 
tice is inhibited outside the trap center due to the large site- 
to- site energy shift. These states were predicted theoretically 
in Ref. [16] and experimentally observed in Ref. lIlTh . Here, 
we show the first dynamic population and in-situ detection of 
such states. Unlike previous experiments, our method allows 
for the transfer of atoms into states that are separated by large 
distances (~ 100 sites) from the original cloud. Furthermore, 
our transport mechanism does not simultaneously broaden the 
spatial distribution unlike tunneling based long distance trans- 
port 1 22]. It also allows for the coherent transfer of a control- 
lable fraction of the cloud in contrast to conventional transport 



using accelerated optical lattices 11231] . 

The condition for the population of these stationary states 
can be understood by returning to Fig.[T] For a given modula- 
tion frequency i^mod, the oscillating wave packets after some 
evolution time have a spatial overlap with localized states at an 
energy difference of hv^nod- At this time a brief modulation 
pulse will induce coupling to these states in correspondence 
with the Franck-Condon principle for transitions between vi- 
brational states. 

Experimentally, we investigate the transfer to localized 
states by adding a second 500 /iS amplitude modulation pulse 
after variable evolution times. The resulting population is 
shown in Fig. [3] a). We have verified that the states are in- 
deed stationary, by checking that their population is indepen- 
dent of an additional hold time in the lattice before in-situ 
imaging. We observe a double peaked resonance structure 
corresponding to the first passage and subsequent return of 
the wave packet following the reflection. In addition, we ob- 
serve a second double peaked resonance after approx. 10 ms, 
corresponding to the wave packet that has completed 3/4 of 
an oscillation after originally traveling downwards. Figure [3] 
shows that roughly 2% of the total population is observed in 
each localized peak, consistent with a transfer of the majority 
of the wave packet population. A determination of the precise 
transfer efficiency and the optimization of the process will be 
a topic of future investigation. 

Figure [3] b) shows the numerical simulation of the pump- 
probe dynamics obtained by adding a second modulation 
pulse to the simulation discussed above. We extract the bound 
state population and plot only the contribution from the out- 
lying positions. The correspondence with the experimental 
data for early times is remarkable, which clearly confirms our 
detailed understanding of the process. Since the numerical 
simulation is performed at 28 kHz instead of 32 kHz the os- 
cillation is slightly faster, which becomes evident as a clear 
shift of the second double peak. 

Finally, the position of the localized states was systemati- 
cally investigated by choosing a modulation duration of 10 ms 
for various modulation frequencies below and above the tran- 
sition frequency z^02 (see Fig. |4l). We have verified that the 
duration is long enough to reach a steady state population. 
The position clearly moves outwards as the frequency is in- 
creased and as the resonance is crossed a new set of close 
lying components appear. The inset of Fig. |4] shows an image 
recorded at a modulation frequency of 36 kHz, displaying the 
main cloud and an inner and outer set of localized states. 

This behavior can, again, be understood based on the spec- 
trum shown in Fig.[Tl At modulation frequencies Vmod ^ ^02 
population can be coupled both to the first and the second 
band, thus resulting in the two observed sets of populated out- 
lying states. The position of each of these states, x^, can be 
obtained by considering the energy balance between the mod- 
ulation, the band, and the harmonic potential energies 
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FIG. 5. (Color online) Individual addressing of localized states, a) 
Population of localized states after a modulation for 10 ms at 30 kHz. 
RF-removal of b) the lowest and c) the main component. 



FIG. 4. (Color online) Positions of the populated localized states as 
a function of lattice modulation frequency. The fits correspond to the 
expected square root behavior (see Eq. ([T]))- 



where i^oo = 0- From this we obtain the expression for the 
localized state position Xi = \/2h{iyrnod — ^oi)/^^trap- Fit- 
ting the data of Fig.lHto this functional dependence allows us 
to attribute the upper and lower branch to the first and second 
bands respectively. From the fit we obtain UqI = 18.5(2) kHz 
and z^Q2 = 33.06(17) kHz in good agreement with first 
(18.1 to 19.3 kHz) and second (30.6 to 37.2 kHz) Bloch bands 
at 5 = 16. This agreement confirms that the process is well 
understood and demonstrates its applicability for the produc- 
tion of localized clouds in the lattice at a long distance from 
the central component. 

The efficient population of localized lattice states offers av- 
enues for e.g. precision interferometric experiments. This 
requires the ability to manipulate the localized states inde- 
pendently. To demonstrate this ability, we have implemented 
selective removal of components as illustrated in Fig. [5] 
Due to the gravitational potential, the condensate is shifted 
slightly downwards compared to the magnetic field zero. The 
gravitational sag in the weak magnetic confinement Xgag = 
^/^trap = 173 fim is large enough to ensure that each cloud is 
located at unique magnetic fields. Selected peaks can thus be 
removed by applying the appropriate RF-field to couple them 
to the unbound magnetic states. This represents an ideal start- 
ing position for future investigations of the coherence proper- 
ties of the localized states. It also enables a new set of experi- 
ments that probe non-linear matter wave dynamics, including 
the reflection and transmission on a BEG 11241] . the collision of 
wave packets, and the investigation of their dispersive proper- 
ties. 

In conclusion, the experimental realization of quasi-free 
wave packets of ultra-cold atoms based on excitations in a 
combined optical lattice and harmonic trap was demonstrated. 
Furthermore, a probe pulse in direct analogy with pump-probe 
spectroscopy allowed for the resonant de-excitation of the 
wave packets into localized lattice states at a controllable dis- 
tance of more the 100 lattice sites from the main component. 
Finally, the independent manipulation of the population of the 
localized states was implemented. 



A particularly promising experimental extension will be the 
use of multiple modulation frequencies. On the one hand this 
will allow for a full analysis of the Franck-Condon type over- 
lap between bound and quasi-continuous states. On the other 
hand it will enable the controlled population of many, spatially 
separated lattice states. It will also be interesting whether the 
transport can be extended to 2D and 3D lattices and com- 
bined with single site addressing ll25n provide the basis of a 
new model of scalable quantum computation. In addition, our 
preparation of localized states in selected bands may serve as 
a particularly pure source for the important investigation of 
the rich structure of extended Bose-Hubbard Hamiltonians in- 
cluding e.g. the supersolid state [26]. The precise control of 
the external degrees of freedom thus makes these wave pack- 
ets a versatile tool both for fundamental investigations and ap- 
plications in the field of quantum technology. 
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